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We investigate novel Landau level structures of semi-metals with nodal ring dispersions. When 
the magnetic field is applied parallel to the plane in which the ring lies, there exist almost non- 
dispersive Landau levels at the Fermi level (Ef = 0) as a function of the momentum along the 
field direction inside the ring. We show that the Landau levels at each momentum along the field 
direction can be described by the Hamiltonian for the graphene bilayer with fictitious inter-layer 
couplings under a tilted magnetic field. Near the center of the ring where the inter-layer coupling is 
negligible, we have Dirac Landau levels which explain the appearance of the zero modes. Although 
the inter-layer hopping amplitudes become finite at higher momenta, the splitting of zero modes is 
exponentially small and they remain almost flat due to the finite artificial in-plane component of the 
magnetic field. The emergence of the density of states peak at the Fermi level would be a hallmark 
of the ring dispersion. 

PACS numbers: 


Introduction .— Semi-metals, usually the reflection of 
unconventional electronic structures at the Fermi surface 
(FS), are related to various anomalous properties and/or 
exotic phases such as unconventional quantum Hall ef¬ 
fect (QHE) in graphene systems [1-4], pressure induced 
anomalous Hall effect in the Weyl semi-metal (SM) [5] 
and non-Fermi liquid phase and peculiar quantum oscil¬ 
lations in the quadratic band-touching SM [6-10]. Also, 
in many cases, they are classified as topologically non¬ 
trivial metals involving surface states [12-15, 18] which 
are generalizations of the concept of the topological in¬ 
sulator [16] to the metallic systems [17]. 

Recently, there have been many suggestions for the 
novel semi-metals with nodal ring FS with different topo¬ 
logical classification schemes [19-27]. Their topological 
non-triviality ensures the existence of surface modes pro¬ 
tected by inversion, time-reversal or certain lattice sym¬ 
metries. Since those candidate materials for the nodal 
ring semi-metal (NRS) are proposed very recently, the in¬ 
vestigations of their physical properties and experimental 
observations are still outstanding open problems. 

In this work, we demonstrate that the NRSs exhibit 
unusual 3D Landau level structures when the magnetic 
field is applied parallel to the plane of the ring. Notic¬ 
ing that the low energy Hamiltonians for various NRSs 
have the same generic structure, we employ the contin¬ 
uum model for Srlr0 3 as an example of the NRS. It is 
explained later that our results are generic and can be 
applied to other materials as well. We show that NRS’s 
Landau levels simulate the adiabatic transition from two 
decoupled graphenes to a Bernal stacked graphene bi¬ 
layer under the magnetic field with an artificial parallel 
component as a function of the conserved momentum. 
During that process, almost flat Dirac zero modes are 


found inside the nodal ring. For some parameters, not 
far from the realistic ones, the 3D quantum Hall effect 
(QHE) may occur with a little doping. Also, we suggest 
that the nodal ring can be probed by the measurements 
of density of states (DOS) under the magnetic field. 

The model for nodal ring semi-metal .— We consider 
the continuum limit of the tight binding (TB) model for 
Srlr 03 near the U point (ky = (0, — tt, n)) in the Brillouin 
zone (Fig. 1 (a)) [28].: 

TL V = 2 t 0 q b T z - tiq b a z T y - t 0 q c v x + f -^-q c (a x - a y )u y T z 
{(hda + t 4 q b )(J x + {t 3 q a + Uqb)tJ v } v z r y 

+t 5 (a x - <Jy)V X T y (1) 

where q = k — ku and a, 6, c represent orthorhombic di¬ 
rections of the lattice. Here, <r a , r Q and v a are Pauli 
matrices where o a is for the J e s = 1/2 Kramers doublet 
and T a (y a ) is for the sublattices B and R ( Y and G)[20]. 
The realistic TB parameters are known as to = —0.6, 
fi = -0.15, t 2 = 0.13, t 3 = -0.2, t 4 = 0.4 and t 5 = 0.06 
in eV. We consider, however, a wide range of the TB pa¬ 
rameters since we are interested in the generic properties 
of NRSs 

TL 11 has both the time-reversal and chiral symmetry 
(C = (JzVyTz). As a result, we have particle-hole sym¬ 
metric doubly-degenerate dispersion relations as follows. 

££c'(q) = C[(v a q a ) 2 + {v b qb) 2 + {v c q c ) 2 + 2tg (2) 

+C' (8 tl{{v a q a ) 2 + ( v c q c ) 2 ) + v A d qlq 2 c ) 2 j 2 

w here v a = R 3 |/\/2, v b = ^(t\ + 2(4f 2 + tf))/2, v c = 
\J (2 tl + tl)/ 2, v d = {4 [v b v c ) 2 - 2tg(t 4 - 2f 2 ) 2 } 1/4 and 
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FIG. 1: (Color online) (a) Low energy band structures on the 
kb = — 7 T plane where the nodal ring (black ellipse) resides, 
(b) The nodal chain appears when 2t2 = 1 4 . Two coupled 
doubly degenerate Dirac cones reside at the points where the 
nodal lines intersect the plane for given q c . In (c) and (d), 
we compare energy spectra around U point calculated from 
the full tight binding Hamiltonian (solid) and the continuum 
model H v (dashed). 


(, (' = ± 1 . Near the U point, they are in good agree¬ 
ment with those of the original TB model as compared 
in Fig. 1(c) and (d). In the continuum model, the nodal 
ring is just an ellipse satisfying 4t 2 = Z§g 2 + ( 2 i 2 + t|)g 2 
with major and minor radii given by r a = 2 ^ 5 /^ 31 and 
r c = 2 |f 5 |/( 2 tQ + f^) 1 / 2 . The ring exists unless t 3 = 0 or 
to = t 2 = 0. Since the dispersion around the ring is linear 
along the radial and kb direction while constant along the 
nodal line, DOS is proportional to energy and vanishing 
at the Fermi level. Interestingly, when 2^2 = £ 4 , we have 
extra nodal lines for |g c | > r c in k a = 0 plane charaterized 
by a hyperbolic curve 4f§ = ^fg+f 2 )*? 2 — 2(4fg+tf-|-2f§)g 2 
as depicted in Fig. 1 (b). We call those consecutively 
occurring nodal lines as the nodal chain. While those re¬ 
sults are valid only around the U point, one can observe 
the nodal chain structure in the full TB model when it 
respects the chiral symmetry {H TB ,C} = 0 which is re¬ 
alized when t xy = td = 0 . 

Numerical analysis of Landau level quantization .— 
Now, we consider the Landau level quantization of the 
NRS. We assume that the magnetic field is applied along 
c-direction so that it is parallel to the plane of the ring. 
We neglect the Zeeman splitting (~ 10” 4 I3[T]eV) [29] 
since it is much smaller than the Landau level spacings 
of this system which is order of 10 meV when B = IT 
as will be shown later. Using the quantization scheme 
q a = (a + o))/\/21-q and qb = *(a — o))/\/21-q, the Lan¬ 
dau level wavefunction is expressed in the form T = 


V°° fr 8 ! r y l r G I c Bi r R I r Y 1 r G -H T « where 1 ; 

Z-m=0 ibi rn ) °n ) °n ; Wi Fn J wirai % 

is the simple harmonic oscillator (SHO) eigenfunc¬ 
tion [10]. If we separate the eigenvector into two 
pieces as A„ = [c® 1 ", c R ^, , c R ^} T and B n = 
[c 4 t, c Gt, c W, c GI] T , they satisfy the following coupled 
secular equations. 


£A„ — \Jn + lMyA n _|_i + L + B„ + y/nM^A n _i, (3) 
eB„ = y/n + 1M_B„+i + L_A„ -I- -^/TzlVI^LB^— 1 , (4) 


where 


m - V itm ~ 1 r c V ( _1 ) mt3 ~* t4 r , x 

om-S, -iwSi- r i"A 5 ) 


m —1 


>Ib 

Jt 2 q c 


n —3 


2 v / 2 /e 


Lf — —toq c — ( 0 (r 45 - r 35 ) - f 5 (r i4 - r^). (6) 


Here, we adopt the representations of the Dirac gamma 
matrices in Ref. [30]. The coefficient is assumed to be 
zero when its subscript n is negative. 

By solving the above equations numerically, we plot 
the Landau level spectra as functions of the conserved 
momentum q c in Fig. 2(a) and (b). Strictly speaking, 
each band is only doubly degenerate, but one can see 
almost four-fold degeneracies inside the ring away from 
the vicinity of its edge at q c = r c . We label them by 
nonzero integer N in such a way that levels with posi- 
tive(negative) energies are marked by positive(negative) 
integers in increasing(decreasing) order from central ones 
near zero energy. 

We observe doubly degenerate two bands which are 
almost flat near zero energy inside the nodal ring ( q c < 
r c ). Those flat bands are fitted nicely by a formula e ~ 
±(1 — ag 2 )“gf exp (bqf) near the ring’s edge as shown in 
Fig. 2(c). The exponent’s coefficient b is found to be 
approximated as 1.886 x as plotted in Fig. 2(d). As a 
result, the energies of central Landau levels (N = ±1,2) 
rapidly reduce to zero as we go inside the ring from q c = 
r c by an amount of 5q c ~ l/(r c Z|). While the splitting 
between those flat modes becomes finite outside the ring 
( q c > r c ), the only exception is when 2t 2 = td for which 
we have four zero modes for any value of q c . 

In addition, we find that, when q c is well inside the 
ring, our Landau level spectra are Dirac-like as shown 
in Fig. 2 (e) and (f). The energies are proportional to 
the square root of the Landau level index and magnetic 
field. On the other hand, the dispersions cannot be fitted 
by the Dirac Landau levels and show linear behaviors 
for large Landau level indices and magnetic field if q c is 
close to the ring’s boundary or outside the ring. However, 
when the nodal chain appears for 2t 2 = td, one can have 
Dirac Landau levels again for larger q c . 

Interpretation via Graphene bilayer under a tilted mag¬ 
netic field — To understand the nature of the peculiar 
Landau level structures of the NRS, we introduce another 
Hamiltonian which is in a block-diagonalized form of 
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FIG. 2: (Color online) In (a) and (b), we plot Landau levels 
of the nodal ring semi-metal as functions of q c for different 
tight binding parameters when the magnetic field (B=1T) 
is along c direction. Blue solid lines are numerical results 
while red dashed curves are analytic ones obtained from 
the fictitious graphene bilayer model near the ring’s cen¬ 
ter. Central Landau levels (N = 1,2) near the edge of the 
nodal ring are drawn by markers (*) and fitted by a formula 
e ~ ±(1 — aq 2 ) a q^ exp (bq 2 ) for various magnetic fields in (c). 
Here, eo = leV is introduced to make the argument of the 
logarithm dimensionless. We show in (d) that the coefficient 
b in the exponent of the fitting function depends linearly on 
-B -1 or l In bottom panels, we show how the energy de¬ 
pends on the Landau level index and the magnetic field for 
given q c . Deep inside the ring ( q c =0, 0.06), we have VNB de¬ 
pendence of the Dirac particle. Those Dirac-like features are 
lost near the ring’s boundary or outside the ring as shown for 
q c = 0.138 and 0.18. In the case of the nodal chain, however, 
we recover the VNB behavior for larger momenta. 


two 4x4 submatrices. Since it gives us exactly the same 
band structures and Landau levels, it can be regarded 
as a unitary transform of H u . Those two submatrices of 
"H d are given by 


n,D 

n K{K') 


V h K(K')) ’ 


( 7 ) 


where 


h /C(K') = V * o.{<la - (Pa)Px ± VbqbPy, (8) 

gg = v c q c (n x - icosOfiy - ism0p z ) (9) 


for q c < r c and 

= v aqaPx ± v b (q b - (p b )p y + (m 0 q c p Z: (10) 
g e = V2t 5 p x (11) 

for q c > r c . Here, p a is the Pauli matrix and the 
plus-minus sign denotes K(K') valley, hj^,^ is the 
Dirac Hamiltonian with the dispersions centered at q = 
(C Pa,0,q c ) for Q c < r c and q = ( 0, (p b ,g c ) for q c > r c 
where p a = ^(2- ( v c q c ) 2 )/vl and p b = vjq c /2v%. 
When q c > r c , we have the mass term with mo = 

— v^/4v 2 = u c sin# while it is massless inside the 
ring. In the mixing term g g, cos 2 9 = v\/A{v b v c ) 2 . The 
transformed Hamiltonian gives us exactly the same en¬ 
ergy spectra in Eq. (3). 

Since h^ K and h^-, are the Dirac Hamiltonians at differ¬ 
ent valleys K and K', one can interpret "H D as a Hamil¬ 
tonian of graphene bilayer with fictitious interlayer cou¬ 
pling g g. Here, £ = + and — correspond to the upper 
and lower layer of the artificial graphene bilayer and we 
denote its basis as i/j = [ca+, cb + , ca_ , Cb_] t . Further¬ 
more, the position of Dirac points in the upper and lower 
layer are shifted in opposite direction. For instance, in¬ 
side the ring {q c < r c ), two Dirac cones are placed at 
q = (pa,0,q c ) and (—p a ,0,q c ) in the q c = const plane 
in the upper and lower layer of the graphene bilayer as 
shown in Fig. 1(b). This kind of shift can be realized 
when the magnetic field is applied parallel to graphene 
layers with strength B || = 2hp a /ed where d is the inter¬ 
layer distance [31, 32]. This fictitious parallel magnetic 
field has its maximum value at the ring’s center ( q c = 0) 
and vanishes at the ring’s boundary ( q c = r c ). The sys¬ 
tem transforms from two copies of graphene monolayers 
to a graphene bilayer with interlayer coupling gg as in¬ 
creasing q c since g g is proportional to q c . 

Applying the magnetic field along the c direction, we 
obtain the Landau level dispersions numerically and find 
that they are exactly the same as the ones from (3) and 
(4). Now, we discuss detailed properties of the Landau 
level spectra of the NRS. 

(i) Deep inside the ring ( q c -C r c ), the Landau levels of 
the NRS at each q c can be considered as those of the four 
independent anisotropic Dirac particles. When q c = 0, 
g g vanishes and we have four decoupled anisotropic Dirac 
Hamiltonians. In this case, the Landau level spectrum is 
e ^ = ±rjlgV2m with four-fold degeneracy where m 
is an integer and vf = ( v a v b ) 1 ^ 2 . To be precise, one 
has exact zero modes only at q c = 0. However, the 
Dirac-like feature is maintained up to quite large mo¬ 
menta where g g is finite because the mixing between 
wavefunctions on two Dirac cones is exponentially small 

(~ e _p “* B ) with the distance between two Dirac cones, 
2 p a , being maximized at q c = 0. In this regime, the 
main role of the finite coupling g g is the renormaliza¬ 
tion of the Fermi velocity. One can project out the 
contribution of the lower layer and construct an effec- 
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tive low energy Hamiltonian around the Dirac point of 
the upper layer by using the resolvent (e — h ^^) -1 
[4, 33]. This gives us a generalized eigenvalue prob¬ 
lem of the form e(I + ( 2v a p 0 )~ 2 )g6gl)'ip + = (h + - 

(2v a po)~ 2 geh.f c ^ K ,^gg)ip + where ip- (_ is the wavefunction 
for the upper layer. For an intuitive picture, let us focus 
on the case 9 <C 1. Notice that the realistic TB param¬ 
eters (t 2 = 0.13, £3 = —0.2 and t± = 0.4) correspond 
to this limit (cos (9 = 0.9969). After a transformation 
to an orthonormal basis set [34], we arrive at a quite 
simple form of the effective Hamiltonian for the upper 
layer, (1 + (v c q c /v a p a ) 2 )~ 1/2 with the renormal¬ 
ized Fermi velocity v'f = (l + (v c q c /v a p a ) 2 )~ 1 ' 2 v f. In 
Fig. 2(a) and (b), we show by red dashed lines that the 
Dirac Landau levels with the renormalized Fermi velocity 
have good agreements with the numerics when g c < r c , 

(ii) We have, as shown in Fig. 2(b), four exact 
zero modes for arbitrary q c when 9 = 0 (2t 2 = 1 4) 
where has only the inter-layer coupling 2 v c q c between 
A_ and B + sites. In this case, our fictitious model is 
exactly the same as the graphene bilayer with Bernal 
stacking under a parallel magnetic field. Although we 
must rely on the numerics to analyze all Landau lev¬ 
els due to the fictitious in-plane component of the mag¬ 
netic field, one can show that we have exact zero modes 
for arbitrary momentum q c . Here, to avoid unneces¬ 
sary complexity, we only provide analytic solutions for 
the case q c < r c [28]. We define the ladder opera¬ 
tor for the upper (£ = +1) and lower (£ = —1) layer 
as = lB(v a (q a - c Pa) + ivbqb) / V^vf which satisfy 
a+ = — \/2valBPa/vf. They have their own SHO 

eigenfunctions u~ and uf which are shifted spatially 
from each other in 6-direction. In this case, one can 
find four zero energy eigenvectors as ip° K x = [u^, 0, 0,0] T 
and ip° K 2 = c 0 [uo ,0,7oUo ,0] T for K valley and ip^, x = 
[0, 0, 0, Uq ] t and ip^, 2 = Co[0, — 70^0,0, Uq ] t for K' val¬ 
ley, where 70 = p a v a /{v c q c ) and c 0 = (1 + 7 q ) -1 / 2 is the 
normalization factor. One can also find analytic form of 
the eigenfunctions at the zero energy for q c > r c in the 
similar way. The reason why we still have zero modes out¬ 
side the ring is that we have four massless Dirac cones 
along the hyperbolic nodal line in k a = 0 plane (Fig. 
1(b)) since mo = 0 for 9 = 0 [35]. 

(to) On the other hand, the four-fold symmetry of 
the zero modes is broken for nonzero 9 although the 
splitting between them is almost negligible. By us¬ 
ing the zero energy solutions inside the ring, we es¬ 
timate their splitting for finite 9 near q c = 0 as 
A « — ^K(ifq(0)]I^Ss:(i(r'), 2 )l = 

d 0 q c ( 1 - aoq 2 ) 1 / 2 exp(-b 0 q 2 ) where a 0 = v 2 /(2tf), b 0 = 
v c^b / v a and do = v c exp(— 2t 2 lft/v 2 ). The splitting of 
the Lowest Landau level is extremely small near q c = 0 
and looks almost flat due to the factor 1^ in the expo¬ 
nent of do which reflects the effect of the huge fictitious 


parallel magnetic field near q c = 0. The derivation of the 
above is only valid around q c = 0 where the mixing from 
higher Landau levels is minimal so that the projection 
onto the central Landau levels is safe. However, inspired 
by this formula, we could obtain the fitting function for 
the bands near the ring’s edge as shown in Fig. 2(c). 

(iv) At q c = r c and 9 = 0, where B\\ = 0, the sit¬ 
uation becomes the usual Bernal stacked graphene bi¬ 
layer under the perpendicular magnetic held and the 
Landau levels are evaluated as = ( , 2 1 ^ 2 v c q c (Xi in + 
((\ 2 ,n + w 4 ) 1 / 2 ) 1 / 2 where A m ,n = 1 + m{2n — l)w 2 and 
u> 2 = v 2 /(2v 2 q 2 lft). The Lanadu levels near the Fermi 


level are described by In the low energy regime, it 
is approximated to « ^v 2 /(v c q c l^)^n(n — 1) which 
is the well-known Landau levels of graphene bilayer [4]. 

3D quantum Hall effect — One of the interesting re¬ 
sults for 9 = 0 case (Fig. 2(b)) is that the zero en¬ 
ergy flat modes are separated from other bands with a 
finite gap. In this case, one can have the 3D QHE by a 
slight doping as indicated by Halperin [36]. He showed 
that the conductivity tensor should be in the form of 
<Jij = e 2 /{2 tt h)etjkGk where eijk is the Levi-Civita sym¬ 
bol and G is the reciprocal vector [36, 37]. When q c = 0 
and q c = r c , the Hall conductance in the ab plane is given 
by alf = 2e 2 /h since our model is equivalent to the 
decoupled spinless graphene layers and Bernal stacked 
graphene bilayer respectively. As we do not have any gap 
closing, the Hall conductances at other momenta are also 
a aF = 2e 2 /6 due to the adiabatic continuation. Then, 
the three dimensional Hall conductance is evaluated as 


_ [ dq c 2 d _ 
ab ~ / 27 t aab ~ 2irh a. 


47T 


( 12 ) 


where a c is the lattice constant along the c direction. 
In the real material, it is expected that 9 0 and also 

there may be other electron or hole pockets at the Fermi 
level. In this case, we expect that the Hall conductance 
would be order of ~ e 2 /ha c even though it is not strictly 
quantized. Further, we expect the strain effect may be 
used to tune the TB parameters close to the ideal case 
of the above. 

Discussion — Although we used a continuum model 
for SrIrC>3, the flat 3D Landau levels at the Fermi en¬ 
ergy can be found in any semi-metals with nodal ring 
dispersion, where the energy spectra are linear along the 
perpendicular directions of the nodal line so that one 
can find Dirac cones for given momentum parallel to the 
ring’s plane as shown in Fig. 1(b). Recently, other can¬ 
didate materials for the NRSs with the above properties 
have been suggested, such as CuaNZn, CuaNPd [21, 26], 
Ca3P2 [22], LaN [23] and so on [24-26]. For example, 
the low energy Hamiltonian of CusNZn, for a given q y 
(q_L = {q x , Qy)), is given by H ~ 2b±q Q (q x - qo)T z +vq r T y 
where q 0 = ±(—Ae/6j_ — q 2 ) 1 / 2 are the positions of two 
Dirac cones. It is also noticed that Ca3P2 has similar 
electronic structures from the band crossings along MYI\ 
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and the linear DOS around the nodal ring. Among these, 
CaaP2 may be the most promising since its ring is free 
from other electron or hole pockets and has a sizable ra¬ 
dius. Flat Dirac Landau levels would appear in these sys¬ 
tems when the magnetic field is applied in the direction 
of the ring’s plane. Because the flat bands have promi¬ 
nent peak in the DOS, one might identify the existence 
of the nodal ring by the scanning tunneling microscopy 
even if it is buried in other dispersive bands. 

This work was supported by the NSERC of Canada, 
the CIFAR, and the Center for Quantum Materials at 
the University of Toronto. We thank Y. Chen, H.- Y. 
Kee, H.-S. Kim and J. Kim for useful discussions. 
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